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1 Abstract 

£f~) ' Effect of Gaussian random magnetic field distribution which is centered at zero on the phase 
t-H , transition properties of isotropic quantum Heisenberg model has been investigated on two (2D) 
and three dimensional (3D) lattices within the framework of effective field theory (EFT) for a two 
spin cluster (which is abbreviated as EFT-2). Beside the phase diagrams and the evolution of 
the magnetization versus temperature curves with the Gaussian magnetic field distribution width, 
critical Gaussian distribution width values, which make the critical temperature zero, have been 
obtained for several lattices. Moreover, it has been concluded that all critical temperatures are of 
the second order and reentrant behavior does not exist in the phase diagrams. 

Keywords: Quantum isotropic Heisenberg model; Random magnetic field; Gaussian 
magnetic distribution 

2 Introduction 

Recently, there has been growing theoretical interest in the random field lattice spin models. The 
model was introduced for the first time by Larkin [I] for superconductors and later generalized by 
Imry and Ma [2]. Ising model which is the most basic lattice spin model with quenched random 
field (RFIM) has been studied over three decades, since this model can be used to describe a wide 
variety of disordered systems. Diluted antifcrromagncts (such as Fe x Zni- x F 2 , Rb 2 Co x Mg 1 ^ x F4 
and Co x Zn\- x F-2) in a homogenous magnetic field behave like ferromagnetic systems in the pres- 
ence of random fields [3] [4]. Structural phase transitions in random alloys, commensurate charge- 
density-wave systems with impurity pinning, binary fluid mixtures in random porous media, and 
the melting of intercalates in layered compounds, such as TiS 2 [5] are the examples of the ex- 
perimentally accessible disordered systems which can be described by RFIM. Besides, RFIM can 
mimic the phase transitions and interfaces in random media [H] [7] , e.g pre- wetting transition on a 
disordered substrate can be mapped onto a 2D RFIM problem j8j . RFIM has also been applied to 
describe critical surface behavior of amorphous semi-infinite systems [9] [10] . 

RFIM has been widely studied in the literature with discrete [TTJ [TU EH HH HH] , as well as 
continous [TBI [TT1 HH HU [20j [21] distributions. Random distribution of the magnetic field may 
produce drastic effects on the phase diagrams and related magnetic properties of the system. 
It has been shown that Ising systems under the influence of discrete symmetric distributions, 
like bimodal lllj and trimodal }12j distributions, exhibit tricritical behavior, while continuous 
symmetric distributions like Gaussian distribution [16] lead to only second order transitions. 

Although the results of the RFIM have been reported for both discrete and continuous distri- 
butions, there has been less attention paid on the random field effects on the Heisenberg model. 
Since Heisenberg model is more realistic model than the Ising model for the spin systems, it is 
important to investigate this model in the presence of quenched random fields. Albuquerque and 
Arruda [22] studied the effect of the bimodal random field distribution on the phase transition 
characteristics of the spin- 1/2 isotropic classical Heisenberg model and they found tricritical be- 
havior within the EFT-2 formulation. Oubelkaccm et al., studied the same system with another 
approach, namely EFT with probability distribution technique and they obtained similar results 
[23] . Albuquerque et al. [24] treated the same system with amorphization effect, again with the 
EFT-2 formulation. Recently, Sousa et al. have studied the effect of the bimodal random field 
distribution on phase transition characteristics of the isotropic classical- and quantum- spin- 1/2 
Heisenberg model within the EFT-2 formulation and also they found a tricritical behavior [25] . 
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All of these works have been restricted to the spin-1/2 isotropic Heisenberg model with bimodal 
random field distribution and they concluded that tricritical behavior exists in this system as in 
Ising model with bimodal random field distribution. They utilized an EFT formalism which is 
characterized by differential operator technique introduced by Honmura and Kaneyoshi for Ising 
systems [26] . 

EFT approximation can provide results that are superior to those obtained within the tra- 
ditional mean field approximation, due to the consideration of self spin correlations which are 
omitted in the mean field approximation. EFT for a typical Ising system starts by constructing 
a finite cluster of spins which represents the system. Callcn-Suzuki spin identities [27| 128] are the 
starting point of the EFT for the one spin clusters. If one expands these identities with differential 
operator technique, multi spin correlations appear, and in order to avoid from the mathematical 
difficulties, these multi spin correlations are often neglected by using a decoupling approximation 
[2"5] . Working with larger finite clusters will give more accurate results. Callen-Suzuki identities 
have been generalized to two spin clusters in Ref. [30] (EFT-2 formulation). This EFT-2 formu- 
lation has been successfully applied to a variety of systems, such as quantum spin-1/2 Heisenberg 
ferromagnet [HI [32] and antiferromagnet [33] systems, classical n- vector model [34] [35], and spin-1 
Heisenberg ferromagnet (36] [37] . 

The aim of this work is to investigate the effect of the Gaussian random field distribution cen- 
tered at zero on the phase transition characteristics of a spin-1/2 isotropic quantum Heisenberg 
model. Quantum Heisenberg model can take into account the quantum fluctuations which dom- 
inates the thermal fluctuations in the low temperatures. Thus, it is expected that it gives more 
reasonable results than the classical one in this low temperature region. We follow the EFT-2 
formulation which is derived in Ref. |31j for this system. 

The paper is organized as follows: In Sec. [31 we briefly present the model and formulation. 
The results and discussions are presented in Sec. [4] and finally Sec. [5] contains our conclusions. 

3 Model and Formulation 

We consider a lattice which consists of TV identical spins (spin-1/2) such that each of the spins has 
z nearest neighbors. The Hamiltonian of the system is given by 



where Si and sf denote the Pauli spin operator and the z component of the Pauli spin operator at a 
site i, respectively. J stands for the exchange interactions between the nearest neighbor spins and 
Hi is the longitudinal magnetic field acting on the site i. The first summation is carried over the 
nearest neighbors of the lattice, while the second one is over all the lattice sites. Magnetic fields 
are distributed on the lattice sites according to a Gaussian distribution function which is given by 



where a is the width of the distribution. According to Eq. ([2]), the average magnetic field value 
of the overall system is zero. The limit a — > covers a pure system, i.e. Heisenberg model under 
zero magnetic field. 

We use the two spin cluster approximation as an EFT formulation, namely EFT-2 formulation 
|31j . In this approximation, we select two spins (namely si and s-i) and treat the interactions 
exactly in this two spin cluster. In order to avoid some mathematical difficulties, we replace the 
perimeter spins of the two spin cluster by Ising spins (axial approximation) |32j . After all, by using 
the differential operator technique with decoupling approximation |29j . we get an expression for 
the magnetization per spin as 




CD 



<i,3> 




(2) 



/// = 



-(si + 4)) = ([A x + mB x ] 



[A y + mB y ] Zo [A xy + mB xy ] Zl ) F (x, y) U=o, y =o 



(3) 
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where s\ and s 2 have zq distinct nearest neighbors and both of them have Z\ common nearest 
neighbors. 

The coefficients are defined by 

A x = cosh ( J Z V X ) B x = sinh ( J Z W X ) 

Ay = COSh ( J z Vy) By = Smh (J Z Vy) (4) 

A X y = COSh [J Z (V a + Vy)] B X y = Sffih [jj, (V X + Vj,)] 

where V 2 - = d/dx and V y = <9/<9y are the usual differential operators in the differential operator 
technique. Differential operators act on an arbitrary function via 

exp (aV x + bV y )G (x, y) = G (x + a,y + b) (5) 

with any constant a and b. The function in Eq. (j3]) is given by 

F{x,y) = jdH 1 dH2P{H 1 )P{H 2 )f{x,y,H 1 ,H 2 ) (6) 

where 

f( tt tt s sinh (f3 X Q ) 

7^,y,xn,jj 2 ; cosh (PX ) + exp (—2/3 J) cosh (f3Y ) [> 

X = (x + y + H 1 +H 2 ), Y =[(2J) 2 + (x-y + H 1 -H 2 ) 2 \ 1 ' 2 , (8) 

with (3 = l/(ksT), ks is the Boltzmann constant and T is the temperature. With the help of the 
Binomial expansion, Eq. (|3|) can be written as 

m = EEE c ;, r mP+,+r ( 9 ) 

p=O g =Or=0 

where the coefficients are 

G' pqr = ( p ) ( g ) ( Z r ) ^-^-^"^W^J/) U=o.«=o ( 10 ) 

and these coefficients can be calculated by using the definitions given in Eqs. (Q| and Let us 
write Eq. ((9]) in more familiar form as 



and 



^C k m k (11) 

and 



zo zo Zl 



Cfe - / ,/ y Sp+q+r,kCp Qr (12) 
p=0<j=0r=0 

where Jjj is the Kronecker delta. It can be shown from the symmetry properties of the function 
defined in Eq. ((6]) and operators defined by Eq. (j4]) that for even k, the coefficient Ck is equal to 
zero. This property is derived in Sec. lAl 

For a given set of Hamiltonian parameters (J), temperature (kgT/ J) and field distribution 
parameter (a), we can determine the coefficients from Eq. (|12p and we can obtain a non linear 
equation from Eq. ([lip . By solving this equation, we can get the magnetization (m) for a given 
set of parameters and temperature. Since the magnetization is close to zero in the vicinity of the 
critical point, we can obtain a linear equation by linearizing the equation given in Eq. (|llj) which 
allows us to determine the critical temperature. Since we have not calculated the free energy in 
this approximation, we can locate only second order transitions from the condition 

d =1, C 3 < 0. (13) 

The tricritical point at which the second and first order transition lines meet can be determined 
from the condition 

d = 1, C 3 = 0. (14) 
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4 Results and Discussion 



The effect of the zero centered Gaussian magnetic field distribution on the isotropic quantum 
Hciscnberg model problem has one parameter as a measure of randomness, namely the width of 
the distribution a. As a increases then randomly distributed magnetic fields with greater absolute 
strengths start to act on the lattice sites. 

In order to see the effect of the random field distribution width (er) on the critical temperature 
of the model, we depict the variation of the critical temperature with a in (ksT c / J, a) plane in 
Fig. Q]for some selected lattices, namely honeycomb (z = 2, Zi = 0), square (z — 3, Z\ = 0) and 
simple cubic (zq = 5, z\ = 0) lattices. As we can see from Fig. ([I]) that increasing er values decrease 
the critical temperature of the whole lattices continuously. This is due to increasing randomness 
in the system: Increasing the randomness via raising the width of the magnetic field distribution 
drags the system to the disordered phase while the spin-spin interaction tends to keep the system 
in an ordered phase. After a critical value of the width, randomness prevails and for the values 
that provide er > cr c , the spin-spin interaction can not maintain an ordered phase even at the zero 
temperature. At a certain value of a, the critical temperature of the simple cubic lattice is greater 
than that of the square lattice, due to the excess number of nearest neighbors of square lattice. The 
same relation is also valid between square and honeycomb lattices. Also due to the same reason, 
simple cubic lattice has a wider ferromagnetic region in a (fee 7c/ J, er) plane than the square lattice 
and the square lattice has a wider one in comparison with the honeycomb lattice. Besides, as we 
can see from Fig. ^ that the whole phase transitions between ordered and disordered phases are 
of the second order and the system does not exhibit a reentrant behavior. 

In Fig. ([TJ (b), we depict the low temperature magnetization versus a curves for the lattices 
mentioned above. Here, we fixed the temperature as fc^T/J = 0.001 and this may be considered 
as the ground state due to the fact that thermal energy supplied by the temperature on the 
system may be neglected in comparison with the spin-spin interaction energy. The behavior at 
this temperature is found to be different from the case in the vicinity of the critical temperatures, 
since increasing a first can not change the ground state magnetization value while it changes the 
critical temperature continuously. After than, the ground state magnetization starts to reduce, and 
vanishes with increasing er. Due to the significant strength of the spin-spin interaction originating 
from the large number of nearest neighbors concerning the simple cubic lattice, the ground state 
magnetization of the simple cubic lattice can become saturated at 1.0 even for large a values, in 
contrast to the square and honeycomb lattices. The a c values at which the critical temperatures 
reduce to zero corresponding to different lattices can be seen in Table (1). As seen in the Table 
(1), the EFT-2 formulation is able to distinguish between the lattices with the same coordination 
number, but different geometry, i.e. the critical width values of the Kagome and square lattices, 
or triangular and simple cubic lattices are different. For instance, er c = 2.262 for Kagome lattice 
is slightly lower than the same value for the square lattice a c = 2.320. Although these two lattices 
have same number of nearest neighbor, this difference comes from the difference between the 
Kagome (zo = 2, z\ = 1) lattice and the square lattice (zo = 3, z± = 0) . 

The critical value of a for the 3D lattices can be compared with the same values for the Ising 
model which are er c = 3.850 for the simple cubic, er c = 5.450 for the body centered cubic and 
a c = 8.601 for the face centered cubic lattices[21]. For the isotropic Heisenberg model these 
critical values are slightly greater than the Ising counterparts. 

Now let us investigate the evolution of the variation of the magnetization curves with tempera- 
ture with increasing randomness in the system. In order to achieve this, we depict the variation of 
the magnetization with temperature curves for some selected values of a in Fig. for honeycomb 
and square lattices. As seen in Fig.JTJ), since the difference between the 2D and 3D lattices regard- 
ing to the behavior of the critical temperature and ground state magnetization is only quantitative, 
we do not expect a qualitative difference between the magnetization profiles of 2D and 3D lattices 
with zero centered Gaussian random field distribution. As seen in Fig. (pj, magnetization versus 
temperature curves related to the honeycomb and square lattices are qualitatively similar. As 
er increases then both critical temperature and the ground state magnetization become reduced. 
However, a = 1.0 curve starts with the same value of the er < 1.0 curves for the square lattice 
while this is not in the case for the honeycomb lattice. This situation supports the behavior of the 
ground state magnetizations given in Fig. ([TJ (b). Moreover, there is not any observed first order 
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Figure 1: (a) Phase diagrams of the isotropic quantum Heiscnbcrg model with Gaussian magnetic 
field distribution in (fcsT c /J, a) plane for selected lattices, (b) Variation of the ground state 
magnetization of the system with (a) for selected lattices. In (b), the temperature has been fixed 
as k B T c /J = 0.001. 

Table 1: The critical value of distribution 
width (cr c ) of the isotropic quantum Hcisen- 
bcrg model with zero centered Gaussian mag- 
netic field distribution. 



Lattice 


zo 


Z\ 


(Tc{H) 


Honeycomb 


2 





1.466 


Kagome 


2 


1 


2.262 


Square 


3 





2.320 


Triangular 


3 


2 


3.875 


Simple cubic 


5 





3.985 


Body centered cubic 


7 





5.563 


Face centered cubic 


7 


4 


8.844 



transition in the curves with increasing a. This also indicates that increasing a can not induce a 
first order transition in the isotropic Heiscnbcrg model. 

5 Conclusion 

In conclusion, the effect of the zero centered Gaussian random magnetic field distribution on the 
phase diagrams of the isotropic quantum Heiscnbcrg model has been investigated in detail. In this 
regard, the effects of the random magnetic fields have been discussed for 2D and 3D lattices. The 
phase diagrams, which are the variation of the critical temperature with the width of the Gaussian 
distribution have been depicted for the honeycomb, square and simple cubic lattices. It has been 
found that there is not any reentrant behavior or first order transition in the system. The critical 
Gaussian distribution width at which the critical temperature of the system vanishes has been 
obtained for several 2D and 3D lattices within the isotropic model. 

Besides, the effects of the increasing randomness of the magnetic field distribution on the 
behavior of the magnetization versus temperature curves have been investigated. We have not 
observed any qualitative difference between the results for 2D and 3D lattices. We have found that 
as a increases then the ferromagnetic region in (m — ksT c / J) plane becomes narrower and finally 
ferromagnetic region disappears right after the critical value of the Gaussian distribution width. 
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Figure 2: (a) Phase diagrams of the isotropic quantum Heisenberg model in (fcgT c /J, a) plane 
with Gaussian magnetic field distribution for selected lattices, (b) Variation of the ground state 
magnetization of the system with (a) for selected lattices. In (b), the temperature has been fixed 
as k B T c /J = 0.001. 

A Symmetry properties of the coefficients 

In order to sec the symmetry properties coefficients given in Eq. (|12j) . let us start with writing 
coefficients given Eq. (flQj) as 

c - "(?)(?)(?) 9 «" (15) 

where 

Q pqr = J dH 1 dH 2 P{H 1 )P{H 2 )@' vqr {H u H 2 ) (16) 



and 



B' pqr (H U H 2 ) = Al^A^A^- r BlBlBl v f (x, y, H U H 2 ) \ x=0>y=0 . (17) 

Here the function defined by Eq. (0 and the distribution function for the magnetic field on the 
sites labeled by 1 and 2 is given by Eq. ([2]). From the definitions given in Eqs. ((4]), ([5]) and with 
using Binomial expansion, ()17[) can be written in the form 



z -p p z -q q Zi—r 



K t . v 



® vqr (H 1 ,H 2 )= ^EEEEE^/K fl ^i,ff2) (is) 

t 1= 0l!i=0t 2 =0D2=0t3=0l>3=0 

where t,v stands for the (t\, t 2 , £3, v\, v 2 , U3) and 

ai = (z - 2ti - 2vi +zi- 2t 3 - 2v 3 ) J z 
a 2 = (z - 2t 2 - 2v 2 +Zi — 2t a - 2v 3 ) J z 

Zo ~P \ ( V \ ( z Q -q \ f q \ f z 1 - r \ / r \ (19) 
ti J\vi)\ t 2 J \v 2 J \ h ) \ v 3 

We can see from Eq. lfl8|) that, expanded form of the expression have both of the terms which 
has / (01, a 2 , Hi, H 2 ) and / (— ai, —a 2 , Hi, H 2 ) for all possible values of ai and a 2 . In other words, 
every term which include / (a±, a 2 , Hi, H 2 ) has a corresponding term / (— a±, —a 2 , Hi, H 2 ). Let us 
focus on these two term. We can see from the definitions given in Eq. (|19|) that the transformation 



(20) 



h 




~P- 


ti 


Vl 


p- 


- Vl 




h 


-> z 


-q- 


t-2 




-> Q- 


- v 2 




h 


-> Zl 


— r — 


h 


V3 


-> r - 


- v 3 
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transforms the terms in Eq. (|19p as 



ai — > — cii 

a 2 -» -a 2 (21) 
K t>v -> (_i)W- 2 (-i+^+^)x tv 

i.e. the terms which have / (ai, a 2 , Hi, H 2 ) and / ( — ai , — a 2 , Hi, H 2 ) in the expanded form of the 
Eq. (|18p have same coefficients if p + q + r is even and same but opposite signed coefficients if 
p + q + r is odd. Thus we can write (fT5)l as 



©pgr (Hi,H 2 ) = 



^^rlf (ai,a 2 ,H 1 ,H 2 ) + f (-ai,-a 2 ,Hi,H 2 )} , p+q+r is even 

t,v 

^2^-[f {ai,a 2 ,Hi,H 2 ) - f (-ai,~a 2 ,Hi,H 2 )] , p+q+r is odd 



t,v 

(22) 

Here, the limits of the sums are identical to those in Eq. (|18p and instead of using six sum as in 
Eq. (|18p. only one sum is used in short notation. 

From Eqs. (0 and (J8]) we can see that the function satisfies 

/ {ai,a 2 , Hi,H 2 ) = -f (-01, -a 2 , -H u -H 2 ) (23) 

then we can write (12^1) as 



©pgr (Hi,H 2 ) = 



X!"^ If (ai,a 2 ,Hi,H 2 ) - f {a x ,a 2 , -H 1: -H 2 )] , p+q+r is even 

t,v 

X!"^ [f {o-i,a 2 ,Hi,H 2 ) + / (a x ,a 2 , -H u -H 2 )\ , p+q+r is odd 

(24) 



K 

t.v 



Thus, we can conclude from Eq. (|24p that, 

& pqr (Hi, H 2 ) = ( " 9 wr ~^ 2) ' P + q t! iSe H e ; • (25) 

P9 ^ S pqr (-Hi,-H 2 ) , p+q+r is odd 

Now, if we look at the integrant of Eq. (|16p . with the help of the Eqs. ([2]) and ([25]), we can 
see that it is symmetric about the origin for the odd valued p + q + r and antisymmetric about 
the origin for the even valued p + q + r in the (Hi, H 2 ) plane. Thus, with using this result in Eqs. 
(fT6f and (fTB"]) then in Eq. (fT2j) we arrive the property 

J = , p+q+r is even - . 

° fe \ ^0 , p+q+r is odd {0> 

and this completes of our derivation. 



7 



References 

[1] A. I. Larkin, Sov. Phys. JETP 31, 784 (1970). 
[2] Y. Imry and S. K. Ma, Phys. Rev. Lett. 35, 1399 (1975). 
[3] S. Fishman and A. Aharony, J. Phys. C 12, L729 (1979). 
[4] J. L. Cardy, Phys. Rev. B 29, 505 (1984). 

[5] Daniel S. Fisher, Geoffrey M. Grinsrcin, Anil Khurana, Physics Today 56, December (1988). 
[6] M.E. Fisher, J. Chcm. Soc. Faraday Trans. 2, 1569 (1986). 

[7] G. Forgacs, R. Lipowsky, T.M. Nieuwenhuizen C. Domb, J. Lebowitz (Eds.), Phase Transitions 
and Critical Phenomena, vol. 14 Academic Press, London (1991), p. 136 

[8] R. Blossey, T. Kinoshita, J. Dupont-Roc, Physica A 248, 247 (1998). 

[9] Y. El Amraoui, A. Khmou J. Magn. Magn. Mater. 218, 182 (2000). 
[10] Y. El Amraoui, A. Hamid, S. Sayouri J. Magn. Magn. Mater. 219, 89 (2000). 
[11] A. Aharony, Phys. Rev. B 18, 3318 (1978). 
[12] D. C. Mattis, Phys. Rev. Lett. 55, 3009 (1985). 
[13] E. F. Sarmento and T. Kancyoshi, Phys. Rev. B 39, 9555 (1989). 
[14] N. G. Fytas, A. Malakis, and K. Eftaxias, J. Stat. Mech. Theory Exp. (2008), 03015. 
[15] I. A. Hadjiagapiou, Physica A 389, 3945 (2010). 
[16] T. Schneider and E. Pytte, Phys. Rev. B 15, 1519 (1977). 
[17] T. Kancyoshi, Physica A 139, 455 (1985). 

[18] Y. Q. Liang, G. Z. Wei, Q. Zhang, Z. H. Xin, and G. L. Song, J. Magn. Magn. Mater. 284, 
47 (2004). 

[19] N. Crokidakis and F. D. Nobrc, J. Phys. Condons. Matter 20, 145211 (2008). 

[20] O. R. Salmon, N. Crokidakis, and F. D. Nobre, J. Phys. Condons. Matter 21, 056005 (2009). 

[21] U. Akmci, Y. Yiiksel, H. Polat, Phys. Rev. E 83, 061103 (2011). 

[22] Douglas F. de Albuquerque, A.S. de Arruda Physica A 316, 13 (2002). 

[23] A. Oubelkacem, K. Htoutou, A. Ainane, M. Saber Chin. J. Phys. 42, 717 (2004). 

[24] Douglas F. de Albuquerque, Sandro L. Alves, A.S. de Arruda Phys. Lett. A 346, 128 (2005). 

[25] J. Ricardo de Sousa, Douglas F. de Albuquerque, Alberto S. de Arruda, Physica A 391, 3361 
(2012). 

[26] R. Honmura, T. Kaneyoshi J. Phys. C 12, 3979 (1979). 

[27] H.B. Callen Phys. Lett. 4, 161 (1963). 

[28] H. Suzuki Phys. Lett 19, 267 (1965). 

[29] T. Kancyoshi, Acta Phys. Pol. A 83, 703 (1993). 

[30] A. Bobak, M. Jascur Phys. Status Solidi B 135, K9 (1986). 

[31] T. Idogaki, N. Uryu Physica A 181, 173 (1992). 



8 



[32] J. Mielnicki, G. Wiatrowski, T. Balcerzak, J. Magn. Magn. Mater 71, 186 (1988). 

[33] Ijanflio G. Araiijo, J. Cabral Neto, J. Ricardo dc Sousa Physica A 260, 150 (1998). 

[34] J. Ricardo de Sousa, Douglas F. de Albuquerque Physica A 236, 419 (1997). 

[35] J. Ricardo de Sousa, Physica A 256, 383 (1998). 

[36] Y. Miyoshi, A. Tamaka, J.W. Tucker, T. Idogaki J. Magn. Magn. Mater. 205, 110 (1999). 

[37] T. Idogaki, A. Tanaka, J.W. Tucker J. Magn. Magn. Mater. 177-181, 157 (1998). 



9 



